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Abstract 

A self-consistent model of the superfluid (SF) state of a Bose liquid with strong inter- 
action between bosons is considered, in which at T — 0, along with a weak single-particle 
Bose-Einstein condensate (BEC), there exists an intensive pair coherent condensate (PCC) 
of bosons, analogous to the Cooper pair condensate of fermions. Such a PCC emerges due 
to an effective attraction between bosons in some regions of momentum space, which results 
from an oscillating sign-changing momentum dependence of the Fourier component V(p) of 
the interaction potentials U(r) with the inflection points in the radial dependence. The col- 
lective many-body effects of renormalization ("screening") of the initial interaction, which are 
described by the bosonic polarization operator II(p,a>), lead to a suppression of the repul- 
sion [V(p) > 0] and an enhancement of the effective attraction [V(p) < 0] in the respective 
domains of nonzero momentum transfer, due to the negative sign of the real part of II(p, ui) 
on the "mass shell" u) = E(p). The ratio of the BEC density no to the total particle density 
n of the Bose liquid is used as a small parameter of the model, no/n <C 1, unlike in the 
Bogolyubov theory of a quasi-ideal Bose gas, in which the small parameter is the ratio of 
the number of supracondensate excitations to the number of particles in an intensive BEC, 
(n — no) /no < 1. A closed system of nonlinear integral equations for the normal En(p, uj) 
and anomalous Ei2(p,tu) self-energy parts is obtained with account for the terms of first order 
in the BEC density. A renormalized perturbation theory is used, which is built on combined 
hydrodynamic (at p — * 0) and field (at p ^ 0) variables with analytic functions £ij(p,e) 
at p — ■> and e — » and a nonzero SF order parameter £12(0, 0) 7^ 0, proportional to the 
density p s of the SF component which is a superposition of the BEC and PCC. In the frame- 
work of the "soft spheres" model with the single fitting parameter — the value of the repulsion 
potential at r = 0, a theoretical quasiparticle spectrum E(p) is obtained, which is in good 
accordance with the experimental spectrum _E cxp (p) of elementary excitations in superfluid 
4 He. 

PACS: 67.57.-z 
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1 Introduction 



An ab initio computation of the spectrum of elementary excitations in the superfluid 
(SF) 4 He Bose liquid remains an actual problem nowadays, despite certain recent 
successes in that direction, like an excellent agreement with experimental data in 
the region of the roton minimum obtained by the Monte Carlo method making use 
of the so-called "shadow wave function" |IJ and by the corellation basic function 
method @] employing modern interatomic potentials for 4 He @-||. At the same 
time, a microscopic field perturbation theory |J-|| calculation of the long- wave 
phonon part of the spectrum E(p) ~ c\p, where c\ is the speed of first (hydrody- 
namic) sound in liquid 4 He, faces principal difficulties. This is due to the fact that 
nonrenormalized perturbation theory gives rise to infrared divergencies and nonan- 
alyticities at p — > and e — > []T2| , which can be cured with the technique of 



combined variables" [|T3j . In the long- wave limit (p — ► 0), those variables reduce 



to the hydrodynamic variables of macroscopic quantum hydrodynamics [13, while 



in the short-wave domain they correspond to the bosonic quasiparticle creation and 
annihilation operators. 

On the other hand, according to numerous precise experimental data on neutron 
inelastic scattering [[HJ-[|n|] and to results in quantum evaporation of 4 He atoms 
fl9|| , the maximal density po of the single-particle Bose- Einstein condensate (BEC) 
in the 4 He Bose liquid even at very low temperatures T < T\ does not exceed 10% of 
the total density p of liquid 4 He, whereas the density of the SF component p s —>■ p 
at T — > pDj . Such a low density of the BEC is implied by strong interaction 
between 4 He atoms and is an indication of the fact that the quantum structure of 
the part of the SF condensate in He II carrying the "excess" density (p s — p ) ^> p 
calls for a more thorough investigation. The questions discussed in this paper are 
both those of the quantum structure of the SF state in a Bose liquid at T = and 
the self-consistent calculation of the spectrum E(p) of elementary excitations in the 
framework of renormalized field perturbation theory f9H- [|l3[ . 



Our approach is based on the microscopic model developed in Refs. pl|— [22|, of 



superfluidity of a Bose liquid with a suppressed BEC and an intensive pair coherent 
condensate (PCC), which can arise from a sufficiently strong effective attraction 
between bosons in some domains of momentum space (see below) and is analogous 
to the Cooper condensate in a Fermi liquid with attraction between fermions near 
the Fermi surface p3| . As a small parameter, one uses the ratio of the BEC density 
to the total Bose liquid density (no/n) <C 1, unlike in the Bogolyubov theory |f24| 



for a quasi-ideal Bose gas, in which the small parameter is the ratio of the number of 
supracondensate excitations to the density of the intensive BEC, (n — n )/n <C 1. 
Because of this, the SF state within the model at hand can be described by a "short" 
self-consistent system of Dyson-Belyaev equations for the normal and anomalous 
Green functions Gik and self-energy parts Eij(p, a;) without account for the diagrams 
of second and higher orders in the BEC density. In this case, the SF component p s 
is a superposition of the "weak" single-particle BEC and an intensive "Cooperlike" 
PCC with coinciding phases (signs) of the corresponding order parameters. The pair 
interaction between bosons was chosen in the form of a finite repulsive potential in 
the "semitransparent" , or "soft" spheres model, whose Fourier component V(p) is an 
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oscillating sign-changing function of momentum transfer p due to mutual quantum 
diffraction of particles. 

As a result of renormalization ("screening") of the initial interaction V(p) due to 
multiparticle collective correlations, which are described by the boson polarization 
operator n(p, cu), the interaction gets suppressed in the domains of momentum space 
where V(p) > 0, and enhanced where V(p) < 0. Such a suppression of repulsion 
and enhancement of attraction is implied by the negative sign of the real part of 
II(p, uj) on the "mass shell" uj = E{j>) for a decayless quasiparticle spectrum. It is 
shown that the integral contribution of the domains of effective attraction in the 
renormalized sign-changing interaction 

v-(p) = v(p) [i - K(p)»n( P , Eip))]- 1 

can be sufficient for the formation of an intensive bosonic PCC in momentum space 
(although not for the formation of bound boson pairs in real space). 

Self-consistent numerical calculations of the boson self-energy, polarization op- 
erator, pair order parameter, and quasiparticle spectrum at T = 0, involving an 
iteration scheme with the single fitting parameter — the value of the repulsion po- 
tential at r = 0, have allowed us to find conditions for the theoretical spectrum 
E{jp) to coincide with the experimentally observed elementary excitation spectrum 
in 4 He |[25||— ||$U||. The roton minimum in the quasiparticle spectrum E(p), which 
corresponds to a maximum in the structural form factor S(q) of a Bose liquid, 
turns out to be directly associated with the first negative minimum of the Fourier 
component of the renormalized potential V(p) of pair interaction between bosons. 



2 Equations for the Green functions and self-energy parts 
in a Bose liquid with a suppressed BEC and an intensive 
PCC in the renormalized perturbation theory 

The main difficulty of the microscopic description of the SF state of a Bose liquid 
with a nonzero BEC is the fact that applying perturbation theory directly || leads, 



as was shown in Refs. P|— 1|12|| , to divergences and non-analyticities at small energies 
e — > and momenta p — > and, as a consequence, to erroneous results in the 
calculations of different physical quantities. Thus, for example, for a Bose system 
with weak interaction, when the ratio of the mean potential energy V(pq)pq (po 
being a typical momentum transfer) to the corresponding kinetic energy p\j1m of 
the bosons is small, the zeroth-approximation polarization operator Il(p, uj) and 
the density-density response function Il(p, uj) calculated to the first order in the 
small parameter of interaction £ = mpo^(^o) <§C 1 ; are logarithmically divergent at 
p — > 0, uj — ► 0, whereas the exact values 11(0, 0) and 11(0, 0) are finite fllQ"| : 

s. dn n n 



where n is the total concentration of bosons, fi the chemical potential, cb = ynVo/m 
the velocity of sound in the Bogolyubov approximation for a weakly nonideal Bose 
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gas PHI , Vq = V(0) the zero Fourier component of the potential, and c the speed of 
sound in the p — ► limit for the spectrum of elementary excitations e(p) ~ c|p| in 
the Belyaev theory ||: 

(2) 



c = V£i 2 (0,0)/m* • 

Here m* is the effective mass of quasiparticles, which is determined by the relation 



1 

m* 



2 
B 



_l_ + d£ n (0,0) <9£ 12 (0,0) 
2m 



d\k\ 



d\k\ 



(3) 



where £n(0, 0) is the normal self-energy part (at k — > 0, e — > 0), and 



B 



de 



£n(0,0) 



<9 2 £ 12 (0,0) Id 



de 2 



+ 2^P^(0,0)r • (4) 



The model, considered in Ref. 0, of a dilute Bose system of hard spheres with a 

small parameter (3 = Jn/ko < 1, in which there is a possibility to exclude the 
infinite repulsion by means of a summation of the "ladder" diagrams, leads to a 
finite value of £12(0, 0) in the zeroth approximation in (3: 



v fn n\ 47ra ° 
£i2(0,0) = n 



m 



(5) 



do being the vacuum scattering amplitude and no the concentration of bosons in 
the BEC (po = mno). 

At the same time, in Ref. taking into account an exact thermodynamic 
equation 



asn(o,o) _ d ni \ =1 _l Si2(0;0) ^o ) 

de \dn J n ' dji ' 



(6) 



where n% = n — no is the concentration of supracondensate bosons, exact asymp- 
totic relations were obtained for the normal and anomalous single-particle Green 
functions: 



Gn(p 0) = -G 12 (p -+ 0) 



n mc 



me' 



c 2 p 2 + id) 



n dn 
m dn 



(7) 



However, it was shown in Refs. P|-||ll| that at p = 0, e = the anomalous 
self-energy part is precisely equal to zero, £12(0,0) = 0. Problems then emerge 
with the determination of the velocity of sound (0) and the asymptotic formulas for 
Gn(p, e) and Gi 2 (p, e) at (p,f)^0| 



C n (p - 0) = -G 12 (p - 0) 



£12(0,0) 



B(e 2 - c 2 p 2 + 18) ' 
because at £i 2 (0, 0) = 0, relations (£D and (|) reduce to identities 

d£n(0,0) 



= 1 



B = 



(8) 



(9) 



so that Eqs. (H) and (H) with account for (||) contain uncertainties of the 0/0 type. 
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With the purpose of fixing these controversies, as well as the infrared divergences 
of Il(p, e) and nonanalyticities in E^-(p, e) at (p, e) — > emerging in the nonrenor- 
malized theory, a renormalization procedure for the field perturbation theory was 
worked out in Ref. [12[], employing the method of "combined variables" fll3| . The 
perturbation theory built on such "adequate" field variables does not suffer from 
infrared divergences at (p, e) — > 0, whose source at T = is the divergence of 
long- wave quantum fluctuations (acoustic Goldstone oscillations). Such oscillations 
are associated with a spontaneous breakdown of continuous gauge and translational 
symmetries in the SF state of a Bose system with a uniform coherent condensate and 
corresponded to the hydrodynamic first sound in liquid 4 He, propagating with the 
velocity of c\ ~ 236 m/s. The choice of combined variables [H| leads to the renor- 
malized anomalous self-energy part £12 (p, e) which does not vanish at (p, e) = 0. 
Then one can formally restore all the results of the nonrenormalized field theory ||- 
0, but now in terms of the renormalized quantities Gik(p,t) and Ejfc(p,e), which 
do not contain singularities at (p, e) — > (save for the pole part Gjfc(p,e) ~ H~ 2 )- 
In particular, the squared velocity of first sound C\ at T — > must be equal to 



Ei2(0,0) 



rrr 



(10) 



where the renormalized effective mass m* is determined by relations @ and (|j) 
with £jfc(0, 0) substituted for Sjfc(0,0) in Eq. (§). In view of the aforesaid, we will 
work with the combined variables [ |l3|j , 



where 



^sh = tp - ^l; 



ty(x 



1 + 



n-L - (n L ) 
2<n L ) 



sh 



^ L (r) 



1 



|k|<feo 



tkr 



IF 



(12) 



(n L )e 



"PL 



Such an approach means that the separation of the Bose system into a macroscopic 
coherent condensate and a gas of supracondensate excitations is made not on the 
statistical level, like in the case of a weakly nonideal Bose gas |?4j , but on the level 
of ab initio field operators, which are used to construct a microscopic theory of the 
Bose liquid. 

The system of Dyson-Belyaev equations which allows one to express the 

normal Gn and anomalous Gu renormalized single-particle boson Green functions 
in terms of the respective self-energy parts E u and S 12 , has the form (Fig. 1): 



Gn(p,e)= G (-p,-e)-E 1 



P, 



/Z(p,e) 



Gi 2 {p, e) = Si 2 (p,e)/Z(p,e) 



Here 



Z(p,e)= G 1 (-p, -e) - E u (-P, -e) G 1 (p, e) - S n (p, e) -|S 12 (p,e 



(13) 
(14) 



(15) 
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Go (P,e) 



P 

e - \- fi + td 

2m 



(6 -> +0) 



(16) 



where /i is the chemical potential of the quasiparticles, which satisfies the Hugengoltz- 



Pines relation [p2 |: 

y u = S 11 (0,0)-S 12 (0,0) . (17) 

Due to a strong hybridization of the single-particle and collective branches of el- 
ementary excitations in the Bose liquid with a finite BEC (no ^ 0), the poles of 
the two-particle and all multiparticle Green functions coincide with the poles of 
the single-particle Green functions G^p, e) [0H§- Therefore the spectrum of all 
elementary excitations with zero spirality is determined by the zeros of the function 
Z(p,e): 



E(p) 



where 



£ + S S U (P,£(P)) - A* 



l 2 



1/2 



|E 12 ( P ,£(p))| 2 + E? 1 ( P ,£(p)), (18) 



£if(p, 



J[E 11 (p,e)±S 11 (-p, 



The (+) and (— ) signs correspond to the symmetric and antisymmetric 
parts of Eu, respectively Relation ensures the acoustic dispersion law for the 
quasiparticles spectrum ( |18"D at p — > with the sound velocity fllCf ). As was shown 
in Ref. f2lfl , for a Bose liquid with strong enough interaction between particles, 
when the BEC is strongly suppressed, one can, when defining Ejfc(p, e) in the form 
of a sequence of irreducible diagrams containing condensate lines, restrict oneself, 
with good precision, to the first (lowest) terms in the expansion over the small BEC 
density (no <C n). Such an approximation is exactly opposite to the Bogolyubov 
approximation for a weakly nonideal Bose gas with an intensive BEC, when 
uq ~ n. As a result, up to terms of first order in the small parameter no/n <C 1, 
for a Bose liquid one gets the "trimmed" system of equations for 
Fig.2): 

E n (p, e) = n A(p, e)V(p, e) + ni V(0) + n(p, e) ; 
E 12 (p, e) = n A(p, e)V(p, e) + # 12 (p, e) , 

where 

d 3 k 



21 



(see 

(19) 
(20) 



%(P 



(27T) 



du 
2^ 



G i3 -(k)V(p -k,e- w)r(p,e,k,w) 



(21) 



\/(p,e) = \/(p)[l-V(p)n( P ,e)]- 1 . (22) 

Here V(p) is the Fourier component of the input potential of pair interaction of 
bosons, V(p,e) is the renormalized ("screened"), due to multiparticle collective 
effects, Fourier component of the retarded (nonlocal) interaction; Il(p, e) is the 
boson polarization operator: 

n(p,6)=i/ 7 ^/^r(p,6,k,a;) 



x {Gn(k,a;)Gii(k- 



J (2vr) 3 J 2tt 
p, e + uj) + Gi 2 (k, a;)Gi 2 (k + p, e + lu)} 



(23) 







T(p, e; k, u) is the vertex part, which describes multiparticle correlations; A(p, e) = 
T(p, e, 0, 0) = T(0, 0, p, e) , and ri\ is the number of supracondensate particles {n\ 
no), which is determined from the condition of conservation of the total number of 
particles: 

. r d 3 k r duj _ . , ni . 
n = n + n\ = n + i / ^-—^ \ — Gn(k, u) . (24) 



(27r) 3 7 2tt 

In the sequel, as well as in Ref. ||21|| , in the integral relations (]21 ) and fl23f ) we will 
only consider the residues at poles of single-particle Green functions Gjj(p, e), ne- 
glecting the contributions of eventual poles of the functions T(p, e, k, uj) and V(p, e), 
which do not coincide with those of Gy(p, e). As a result, taking into account re- 
lations fll3f )-(pOD, Eqs. (|2l|) on the mass shell e = E{p) assume the following form 
(at T = 0): 



*u(p,#(p)) 



d 3 k 



2 J (2?r) 
xV(p-k,E(p 



3 T(p,E{p)-k,E(k)) 



E(k)) 



A(k,E(k)) 
E(k) 



(25) 



*ia(p,£(p)) = - 



d 3 k 



r(p, E(p); k, £?(A;))V'(p - k, E(p) - E(k)) 



where 



E(p) 



+ 



2 J (2?r) 3 

n A(k, E(k))V(k, E(k)) + tt 12 (k, E(k)) 
E(k) 

A 2 (p, E(p)) - [n A(p, E(p))y (p, E{p)) + * 12 (p, £(p)) 
1 



1/2 



*ii(p,£(p))-*ii(-p,-£(p)) 



(26) 



(27) 



A( PjJ E7(p)) = n A(p,E(p))V(p,E(p)) 
1 

+ 2 



*n(p,£(p))+*ii(-p,-£(p)) 



*ii(0,0)+*i 2 (0,0) + ^- . 

2m 

(28) 



The total quasiparticle concentration in the Bose liquid is determined by the relation 



n 



n + 



d 3 k 

J2mf 



A(k,E(k)) 
E(k) 



- 1 



(29) 



From Eqs. ( pf) and fl28|) it follows that the quasiparticle spectrum E(p), because 
of the analyticity of the functions ^(p, e), is acoustic at p — ■> 0, and its structure 
at p depends essentially on the features of the renormalized pair interaction of 
bosons. The theoretical spectrum (B7|) must be close to the experimental spectrum 
of elementary excitations in the 4 He Bose liquid [^5[| -[0O[| if this model is to be 
applicable for the description of the SF state in 4 He (see below). 



7 



Note that when the BEC is totally absent (no = 0), equation ( p5|) becomes 
homogeneous and degenerate with respect to the phase of ^^(p). It is then akin to 
the Bethe-Goldstone integral equation for a pair of particles in momentum space 

*( p) = -/w np " k) ^(IFn- < 30 > 

with zero binding energy, Q = 0, which has a nontrivial solution \&(p) 7^ only 
when there is attraction V(q) < in a broad enough region of momentum transfer 
q. By virtue of this analogy, the function ^12 (p) at no = can be taken to be 
the PCC order parameter [2I|— [22|], which describes condensation of boson pairs in 



momentum space (identical to the Cooper condensate of fermion pairs 

The degeneracy of equation (^) with respect to the phase of \&i 2 (p) at n — > 
allows for the condition of stability of the phonon spectrum c\ = \E'i2(0)/m* > 
to be met by means of choosing the appropriate sign (phase) of the pair order 
parameter ^12(0) > 0. However, the model of the SF state with PCC and no BEC 



21fl-||22|| implies a few paradoxes, such as a finite energy gap in the single-particle 
spectrum at p — 0, exponential asymptotic behavior of the pair correlation function 
^(r)^(r / )^) at |r — r'| — > 00, half-integer quantum of circulation of the SF velocity 
k = h/2m etc. 

Indeed, at p —>■ and n Q 7^ 0, Eq. fl26|) , due to the isotropic momentum de- 
pendence of the_ spectrum E{p) and the functions V(p) = V(p, E(p)), A(p) = 
A(p, E{p)) and ^12 (p) = ^^(Pj E(p)) reduces to the form 



DC 



*12(0) = ~j^y 2 J [n A 2 (k)V\k) + A(k)V(k)* 12 (k)} . (31) 

The first integral addend on the right-hand side of Eq. (|3l|) being always negative, 
the value of ^12(0) can be negative as well. The condition ^12(0) < means that 
the phase of the PCC is opposite to phase of the BEC, because n > 0. Moreover, 
in this case, in spite of the condition A(0)V(0) > (which ensures that the system 
is globally stable against a spontaneous collapse), at sufficiently small densities of 
the BEC, in accordance with Eq. fl20|), the values 



Si 2 (0, 0) = n A(0)V(0) - \^ 12 (0)\ (32) 

become negative if | (0) | > noA(0)V(0), which corresponds to an instability in 
the phonon spectrum (c\ < 0). However, if the pair interaction between bosons in 
a broad enough region of the momentum space has the character of attraction, i.e., 
A(k)V(k) < at k 7^ 0, and if the magnitude of that attraction is large enough 
(see below), the second (positive) addend on the right-hand side of Eq. (|31| ) can 
outweigh the first (negative) one if the BEC density is small enough (no n). 
Then ^12(0) will be positive, and the phase of the PCC will coincide with phase of 
the BEC, so that E 12 > and c? > 0. 

Since at T = the density of the SF component p s , on the one hand, coincides 
with the total density p = ran of the Bose liquid and, on the other hand, p s 
is proportional to £12(0,0) which plays the role of the SF order parameter, with 



S 



account for and fl3~T|), one gets the following relations: 



p s = Po + Ps = /3m 



Si2(0,0) 
A(0)V(0) 



/3m [n (l - 7) + *] 



where 
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(2tt) 2 A(0)F(0) 



fc 2 dfc r 



(2vr) 2 A(0)\/(0) / 



k 2 dk 



A(k)V(k) 



A(k)V(k)* 12 (k) 



(33) 



(34) 



(35) 



1 — 7) 1 . This means that the density 



and (3 is a certain dimensionless constant. Since the density of the single-particle 
BEC is equal to p — rnn , we obtain (3 
of the "Cooperlike" PCC is 



p s = mrii 



m\I/ 



7) 



(36) 



n 



n is then determined from relation ( p9|) , and 



where the concentration n\ 

for liquid 4 He at T — ► 0, in accordance with the experimental data 
should be approximately 90% of the full concentration of 4 He atoms in liquid helium 
n = 2.17 ■ 10 22 cm -3 . Thus, the SF component of Bose liquid in this model at T = 
is an effective coherent condensate [12] which is a superposition of the weak one- 
particle BEC and intensive PCC. 



3 Choice of the pair interaction potential in the Bose liquid 



To describe interaction of helium atoms in real space, various semi-empirical po- 
tentials are conventionally used, which describe strong repulsion at small distances 
and weak van der Waals attraction at large distances. However, most of those po- 
tentials are characterized by a strong divergence at r — > 0, like, for instance, the 
Lennard- Jones potential 

12 



U LJ (r) 



r > r r 



(37) 



Such potentials are not suitable for the description of pair interaction in momentum 
space, since the respective Fourier components 



V(p) = J d 3 r U (r) exp (ipr) = — J rU (r) sin (pr) dr 



(3? 



are infinite, diverging at the lower limit. Lately, in the calculations of interatomic 
interaction and possible bound states, i.e., He2 molecules, one uses more up-to-date 
potentials, like the Aziz potential 0-0: 

Aexp(— ar — Pr 2 ) — F(r,r )(cer~ 6 + Csr~ 8 + Cior~ 10 ) , 



U A (r) = 
F(r, r ) 



exp [-(r /r - l) 2 ] ELo > r < r o 
1 , r > r 



(39) 
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where A = 1.8443101 x 10 5 K, a = 10.43329537 A" 1 , (3 = 2.27965105 A~ 2 , 
c 6 = 1.36745214 K x A 6 , c 8 = 0.42123807 K x A § , c w = 0.17473318 K x A 10 . Such 
potentials remain finite at r = due to the nonanalytic exponential r dependence, 
which suppresses any power divergence at r — > (see Fig. 3). However, employ- 
ing its Fourier component in solving the nonlinear integral equations (f25|)-(|2"6|) is 
technically difficult while not conclusive by itself: In an intrinsically many-body 
problem like the one at hand, collective effects are certain to play an essential part 
and to render the subtleties in the shape of the two-body potential largely irrele- 
vant. In order to retain the crucial features of the system yet not be overwhelmed 
by technical complications, we will utilize a model potential of "soft spheres" de- 
scribing (unlike the "hard spheres" model [3~3~|) finite repulsion in a certain bound 



region, which accounts for effects of mutual quantum diffraction of bosons in the 
Bose liquid. 

In this context, consider a model potential in the form of a Fermi-type function 
in real space (Fig. 4a) 

^ F (r) = ^exp(^-^+lj , (40) 

which at b = degenerates into a "step" of finite height Uq at r < a. In this latter 
case the Fourier component is expressed in terms of the first order spherical Bessel 
function (see Fig. 5): 

v(p) = v ^; Jl (x) = sin(x) - 2 XCQS(X) - (4i) 

pa x 2 

where Vq = 3V^(0) = AnUoa 3 . The same oscillating Fourier component is charac- 
teristic of a smooth potential V(r) in the form of a Lindhardt-type function p3| . 
having an infinite negative derivative at the inflection point r = a (see Fig. 4b): 



(1 — r 2 /a 2 ) , a + r 



1 + r — -In 

2r/a 



a — r 



(42) 



Formally, this problem is an inverse to the one of periodic oscillations of spin density 
in real space of interacting spins in metal (so-called Ruderman-Kittel-Kasui-Yoshida 
oscillations f3~4| ). 

Note that the amplitude of oscillations of the Fourier component for the Fermi 
type potentials ([!(]) at b ^ is damping exponentially with the increase of the 
parameter b, due to the decreasing absolute value of the negative derivative at the 
inflection point (see Fig. 5, inset). 

Such oscillations of the Fourier component of the pair potential U (r) in momen- 
tum space arise even in the absence of attraction in real space and are a consequence 
of quantum diffractional effects of mutual scattering of the particles. This means 
that the existence of negative values V(p) < 0, i.e., of effective attraction in some re- 
gions of momentum transfer p, is not directly associated with van der Waals forces, 
which are explicitly taken into account in the sign-changing (with respect to r) 
Lennard- Jones or Aziz potentials (see Fig. 3). 
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If one substitutes the oscillating potential PIQ into the Bogolyubov spectrum of 
a dilute quasi-ideal Bose gas 



+ 2nV(p) 

2m 



1/2 



(43) 



then, by choosing two parameters, V and a, independently, one can achieve a 
rather satisfactory coincidence of the spectrum Eg(p) with the elementary excitation 
spectrum E exp (p) in liquid 4 He derived from neutron scattering experiments p5| - 
f29fl (Fig. 6, solid and dotted curves). However, the spectrum (|43|) with the potential 
(P|) turns out to be unstable at large values of Vo, because Eg(p) < in some range 
of p (Fig. 6, dashed curve). Moreover, the Bogolyubov model of the quasi-ideal Bose 
gas with an intensive BEC at T — > (uq — > n) is not applicable to the description 
of the Bose liquid with a strongly suppresed BEC (n <C n). 

On the other hand, multiparticle collective effects in the Bose liquid, according 
to Eq. fl22|) , lead to an essential renormalization of the pair interaction, which 
determines the normal and anomalous self-energy parts, Eqs. (|i~9|) and (p0|). An 
important feature of the renormalized interaction fl22|) is that in the regions of 
phase volume (p, a;) where the real part of U(p,uj) is negative, the repulsion (when 
V{p) > 0) gets suppressed while the attraction (when V(p) < 0) gets effectively 
enhanced. This fact was first noted in Ref. p5j and used in Ref. |36| where the 
integral equations ( p5|) and (26) were solved with the seed potential of the hard 
spheres model (see (0)). Note that in Ref. PEfl , a possibility for bound pairs of 



helium atoms to form not only in momentum space, but also in real space was 
discussed, which allowed one to interpret the anomalously large effective mass 
of the dope 3 He atoms in the SF Bose liquid 4 He P7|-P5fl as the mass of a bound 
3 He- 4 He pair, equal to M = m 3 + m 4 = 7m 3 /3. 

In this paper, we take into account the explicit momentum dependence of the 
real part of polarization operator Il(p, E(p)), which can be represented as (see 
Appendix): 



x 



»n( P) s(p)) = iy 



d 3 k 



F(p,k) 



(2n) 3 E(k) -E(\k 



-p\)-E(p) 
^+(k, P ) 



(44) 



E(k)[E(k) + E(\k — p\) — E 



E(k)[E(k) + E(\k — p\) + E 



where 



^-(k,p) 



E(k) - E(p) 



E(k) + — -fi + i:n(k,E(k)) 
2m 



(k- P ) 

2m 



-H + t n (k-p,E(k)-E(p)) 



(45) 



+ S 12 (k, E{k))t l2 {k - p, E(k) - E{p)) 
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E(\k - p|) + - /i + E n (k - p, E(\k - p|)) 



£(|k - p|) + E(p) + + Sn(k, E(\k - p|) + £(p)) 



(46) 



+ S 12 (k - p, E(\k - p|))E 12 (k, £(|k - p|) + E{p)) . 

As follows from Eq. (|44]), if the quasiparticle spectrum -E(p) is stable with respect 
to decays into a pair of quasiparticles |j/fj,P5|, i.e., if for all p and k the following 
conditions are fulfilled: 



E{p)<E(k) + E(\k-p\ 



E(k)<E(p) + E(\k-p\), 



(47) 



the common denominator in front of the curly braces is always negative, 

E(k) — E(\k — p|) — E(p) < , (48 



whereas the denominator of the first term in the curly braces is always positive, 



E{k) + E{\k-p\)-E(p) >0 
and smaller than the positive denominator of the second term 

E(k) + E(\k-p\)+E(p) >0. 



(49) 



(50) 



This means that the integrand of Eq. (H) is negative if the functions F±(k, p) are 
positive. According to the numerical calculations ||36|| , [BO] in the framework of the 



"hard" and "soft" sphere models, F±(k, p) > for all k and p, so that the real 
part of n(p, E(p)) is negative, because T > (see below). One should note that the 
actual experimental spectrum of elementary excitations in liquid 4 He is decaying at 
small enough momenta. However, this will not change the negative sign of n(p) 
due to the integral character of expression ( [4*4]) (cf. ^0[). Note also that here we do 
not take into account the imaginary part of n(p, u), which determines the damping 
of quasiparticles and the dynamical structure factor (see Appendix). 



4 The iterative scheme of calculation of the quasiparticle 
spectrum 

In order to calculate the quasiparticle spectrum within the model of a Bose liq- 
uid with a suppressed BEC and intensive PCC being considered, at first, using 
Eqs. (^) and (p6[), a numerical calculation in the first approximation of the func- 
tions $i(p) = \&n(p, -E'o(p)) and ^i(p) = ^i2(p? -^o(p))? was conducted. For the 
zeroth approximation, the Bogolyubov spectrum E (p) = Eb{p) and the "screened" 
potential (1221) with account for Eq. (f4"l"D were taken: 



Vo(p) = . . r (51 

pa - V Uoji(pa) 

at some constant negative value of n . Then, using the functions $i(p) and ^i(p) 
obtained, the first approximation for the polarization operator IIi(p) was calculated, 
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using Eqs. (|P|)-(|4"6l) at T — 1. Here, too, the Bogolyubov spectrum (^), which is 
the best fit to the empiric spectrum E exp (p) for liquid 4 He (see Fig. 6) was chosen 
as the zeroth approximation for E(p). The limiting value (0) was compared with 
the exact thermodynamic value of the polarization operator of the 4 He Bose liquid 
at p = and u = [cf. ([TJ)], which determines the compressibility of the Bose 
system: 11(0,0) = — n/mcf. 

The absolute value |I1(0,0)| turned out to be almost 1.5 times greater than the 
calculated value |IIi(0)|. This provides an estimate of the vertex r at p = in the 
first approximation as r\ = Ax — 1.5. The second approximation $2(p) an d ^(p) 
was obtained from Eqs. (|25"D , (|26"D with the constant value r x = Ai and the first 
approximation for the renormalized potential, Eq. 



KlW = P a-V U l{p M P a) ■ (52) 

Such an iterative procedure was repeated four to six times and used to improve 
precision in the calculation of the polarization operator. At each stage, equations 
(P7|) and fl2"8"|) were used to reproduce the quasiparticle spectrum E(p), and the rate 
of convergence of the iterations was watched, as well as the degree of proximity of 
E{p) to the empirical spectrum E exp {p). 

The only fitting parameter in these calculations was the amplitude Vq of the 
seed potential (|4lD at the value of a = 2.44 A, which is equal to twice the quantum 
radius of the 4 He atom. The BEC density, in accordance with the experimental 
data Hl5f-||l9ll, was fixed at n = 9%n = 1.95 ■ 10 21 cm -3 . The computation has 



resulted in a quite satisfactory agreement of the theoretical spectrum E(p) with 
E exp (p). Figure 8 depicts the momentum dependence of H(p) as obtained with five 
iterations, while Fig. 9 shows the self-consistent p dependences of $(p), ^{p), and 
A{p) obtained from Eqs. (0), (0), and (j2"8"|). One notices that the functions $(p) 
and H(p) are negative at all p, whereas ^(p) and A(p) are positive. The locations 
of the deep minima of $(p) and A(p) practically coincide with the location of the 
minimum of the potential ( f22"| ) (see Fig. 7). 

Finally, in Fig. 10, the solid curve is the theoretical quasiparticle spectrum E(p) 
obtained from Eq. (p7|), and the dots are the experimental spectrum obtained from 
data on inelastic neutron scattering in 4 He p5|-||29||. In the calculation of E(p), 
the fitting parameter Vo was chosen in such a way that the phase velocity of quasi- 
particles E{p)/p at p — > coincide with the speed of the first hydrodynamic sound 
Ci ~ 236 m/s in liquid 4 He. This value of Vo corresponds to the repulsion potential 
of the "soft" spheres model U = Vo/(4vra 3 ) = 1552 K at a = 2.44 A. We see that 
there is a satisfactory agreement of E(p) with E exp (p) in the region p < 2.5 A . 

For p > 2.5 A , the theoretical spectrum E[p) lies somewhat higher than E exp (p), 
which, apparently, has to do with the fact that the vertex function T(k, p), a de- 
creasing function of p, was replaced with a constant value T ~ 1.5 for all p. 

Of course, the value of E(p) at large momenta should not exceed the doubled 
value of the roton gap A r = 8.61 K lest the spectrum becomes decaying B9]. To 



check this, we have approximated the vertex A(p) = T(0,p) with a slowly decreasing 



function, falling down from f = 1.5 to T = 1.1 on the interval 2.1 A < p < 
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3.8 A . The resulting theoretical spectrum is shown in Fig. 11, together with the 



experimental data f25j-[|9j (light circles) and the latest results |3(J (asterisks) of 
measurements of the spectrum at 2 A < p < 3.6 A , beyond the roton minimum. 
Evidently, such an approximation for the vertex part yields a much better agreement 
of the theoretical and experimental spectra at large values of momentum. Note 
also that the self-consistency of this model is corroborated by the fact that the 
theoretical value of total particle concentration calculated from Eq. (^), = 2.12- 
10 22 cm -3 , is quite close to the experimental value for liquid 4 He, n = 2.17-10 22 cm -3 
(at n = 9%n). On the other hand, the concentration rii of supracondensate 
particles, calculated from Eqs. fl34])-(|36"D at the values of the parameters indicated, is 
about 0.93 n, which is also in good accordance with experiment, taking into account 
that the BEC density is determined up to ±0.01 n. Also, when formulating the 
approximate theoretical model, quadratic terms in the small parameter n /n <C 1 
in Eqs. (|ITf)-(pO"l) were omitted, which also introduces an error of the order of 1% 



5 Conclusions 

Thus, the model of the SF state of a Bose liquid with a single-particle BEC sup- 
pressed because of interaction and an intensive PCC, based upon a renormalized 
field perturbation theory with combined variables [II|— [I3|, allows one to obtain a 



self-consistent "trimmed" system of nonlinear integral equations for the self-energy 
parts SjjQo, e), by means of truncating the infinite series in the small density of the 
BEC (n /n <C 1). By the same token, one can work out a self-consistent micro- 
scopic theory of a superfluid Bose liquid and perform an ab initio calculation of 
the spectrum of elementary excitations E(p), starting from realistic models of pair 
interaction potential U(r) possessing finite Fourier components. It is shown that 
for a repulsive potential in the framework of the "soft spheres" model, the Fourier 
component V(p) is an oscillating sign-changing function of momentum transfer p. 
This means that in certain regions of momentum space at p ^ there is an ef- 
fective attraction between bosons, V(p) < 0, which has nothing to do with van 
der Waals forces and has a quantum mechanical diffraction nature. That attraction 
gets substantially enhanced due to multiparticle collective effects of renormalization 
("screening") of the initial interaction, which are described by the boson polariza- 
tion operator II(p, u). The enhancement of the attraction happens because on the 
"mass shell" u = E(p), the real part of n(p, E(p)) is negative in the whole region of 
momentum where the quasiparticle spectrum E(jp) is stable with respect to decay 
p9| . It is necessary to emphasize that this negative sign of 9ft II(p, E{p)) is only 
characteristic of Bose systems, in which the single-particle and collective spectra 
coincide with each other and are measured from the common zero of energy, unlike 
the Fermi systems, in which the single-particle excitation spectrum begins at the 
Fermi energy, due to the Pauli principle. Therefore, in the 3 He Fermi liquid there 
can be no corresponding effective enhancement of the negative values of the same 
"input" interaction potential V(p), so that the formation of Cooper pairs is only 
possible for nonzero orbital momenta, due to the true weak van der Waals attrac- 
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tion between fermions |?I]]-[|42|]. Apparently, it is this fact that has to do with the 
critical temperatures of the SF transition in 4 He and 3 He differing by three orders 
of magnitude. 

The rather strong pair attraction of bosons in momentum space for 3? H(p, u) < 
forms an intensive PCC, which, together with a weak BEC, constitutes a single 
coherent condensate, making up the microscopic foundation of the SF component 
of the Bose liquid p s ~ Ei 2 (0, 0). On the other hand, the oscillating nature of 
the renormalized Fourier component of the potential V(p) (see Fig. 5 and Fig. 7) 
leads to a nonmonotonic behavior of momentum dependences of the mass operators 
Sn(p, E{p)) and Ei 2 (p, and, as a consequence, to the emergence of a roton 

minimum in the quasiparticle spectrum E(p), which is directly connected with the 
deepest first negative minimum of V(p). 

It is necessary to emphasize that the amplitude of the "soft spheres" model re- 
pulsion potential obtained, Uq = 1552 K, which corresponds to very good agreement 
between the theoretical quasiparticle spectrum E(p) and the experimental spectrum 
of elementary excitations in 4 He is smaller than the value of the Aziz type potential 
at r — > (see Fig. 3, inset). It is a result of strong quantum diffraction effects in 
the Bose liquid, because the average distance between particles is of the order of or 
less than the de Broglie wavelength of the bosons. 
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7 Appendix 

The polarization operator ([19]) can be calculated without account for the vertex 
part T, making use of expressions (H)-(H) m the form 



/d k 
■j—^ [L u (p, k, u) + L 12 (p, k, u)\ , (A.l) 

where 



/dz ~ 
—G ij (k,z)G ij (k-p,z-uj). (A.2) 

Assume that the Green functions Gij have only one pole within the integration 
contour and are equal to 

^" ^-i^U ( ^ 0) - (A - 4) 
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Calculating the integrals (A. 2) with account for the poles at the points e = E(k) 
and e = E(\k — p|) + uj yields 



Ln(p,k,a;) 



2[E(k)-E(\k-p\) 
E{k) - ^ 



UJ 



E(k) + — -fi + J:n(KE(k)) 



uj 



2m 



-ti + Vu(k-p,E(k) 



E(k) [E(k) + E(\k — p|) — uj] 



£(|k- P |) + 



2m 



-V + X n {k-p,E(\k-p\)) 



x 



S(|k-p|)+w+i--/i + Ei 1 (k,S(|k-p|) 



E(\k-p\)[E(k) + E(\k-p\)+uj] 




(A.5) 



Li 2 (p,k,o;) 



2 [£(*;) -£(|k-p|) 



uj\ 



£ 12 (k,£(fc))£ 12 (k-p,E(fc)-^) 
[ J B(Jfe) + J B(|k-p|)-w] 



E 12 (k, £(|k - p|) + u;)£ 12 (k - p, £(|k - p|)) 



£(|k-p|)[I<;(A;) + £(|k-p|) 



uj 



In the statistical limit (uj 

Ln(0,k,0) 
9 / 

+ 



0,p 
1 



0), expression (A.5) reduces to 



(A.6) 



e(k) 




1 2 



me(k) M 

ok 



e(k) + — - t x + t 11 (k,e(k)) 
2m 



e(k) + - fi + E u {k,e{k)) 

2m 



(A.7) 



It follows that in a large region of momentum space, /n(0, k, 0) < 0. The same 
result is obtained for the function (A.6) at p = and uj = 0, i.e., Li 2 (0, k, 0) < 0, so 
that the static bosonic polarization operator 11(0, 0) is negative, which corresponds 
to a suppression of the "screened" repulsion at p — > 0. From Eqs. (A.5) and (A.6) 
it can also be seen that on the "mass shell" uj = E(p), the integrals L n and L 12 
remain negative in a wide region of momentum space because of the negative sign 
of the common denominator E(k) — E(\h — p|) — E(p) < and the positive sign of 
the denominator E(k) +E(\ln — p|) —E(p) > due to the fact that the quasiparticle 
spectrum E(p) is decayless [conditions fl47Dl- 

Thus, the real part of the polarization operator (A.l) at uj = E(p) is negative on 
the whole range of p. In order to determine the imaginary part of II(p, uj), which 
describes Landau quantum damping of bosons, one has to calculate the main value 
of the following integral: 



L(p,k,w) = 
{[e + a(k,e)] [e 



i 

2^ 

- uj - 



V.p. / 

— oo 



de 



[e 2 -E 2 (k)] [(e-uj) 2 -E 2 (k-p)} 
p, e - uj)] + S 12 (k, e)E 12 (k - p, e - u)j 



(A. 
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where 



o-(k,e) = — -/i + S n (k,e) 



Factoring the denominators in the integrand of(A.8) into simple fractions, one 
obtains 

oo 

L(p,k,u) = 27r ^ ( |^_ p | ) V -P- / de[M + (p,k,w,e)-M_(p,k,w,e)] (A.9) 



where 



„, / , , x eC±(p,k,u,e) + D ± (p,k,u,e) C±(p,k,o;,e) 
M±(p, k, w, e) = — — 1 ^ ; (A.10) 



e 2 -E 2 (k) 



e-0JTE(\k-p\) 



C±(p, k,w,e) 



{[£(*;) + a(k,e)]a(k- P ,e-cc>) 



£ 2 (A;) - [w ± £(k - k)] 2 
+ [E(\k - p|) + [<r(k - p, e - u) - u\ + t 12 (k, e)E 12 (k - p, e - u)} ; 



(A.11) 



£>±(p, k, u, e) = cr(k, e) + a(k - p, e - u) - uj + [u ± E(\k - p|)] C±(p, k, e, a;) . 

(A.12) 

If one neglects the explicit e dependence of £y(k, e), the integral (A.9) vanishes. 
On the other hand, from Eq. (A. 10) it follows that the nonvanishing contribution 
into L(p, k, a;) is given by the odd in e parts of the functions C±(p, k, u, e) and the 
even in e parts of the functions D±(p,k,u,e). Since Ei 2 (p, k) is an even function 
of e, and £n(p, k) contains both an even Sf 1 (p, k) and an odd E° 2 (p, k) part, the 
expressions for C± and D±, which provide for nonzero values of QT[(k, cu), can be 
cast into the form 



C±(p,k,u,e) 



E 2 (k) -[u± E(k - k)]' 



{^(k-p.c-w) 



X 



E(k) + E(\k - P |) + uj + S s n (k, e) - fi + — 



(A.13) 



+ £n(k,e) 



Sn(k-p,e-cj) + 



(k- P 

2m 



,21 



At(p, k, w , e) = S s n (k, e) + S s n (k - p, e - - 2// 



2m 



+ [w±£?(|k-p|)] \t s u (k-p,e-u;) 



(k- P ) ; 

2m 



— 



+ £(|k - p|) + to + £^(k, e) - n + 



2m 



+ t a n (k, e)t a n (k -p, e -u) + S 12 (k, e)E 12 (k - p, e - w)} . 



(A.14) 
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At the same time, a nonzero value of S n(p, u) implies that the retarded renormal- 
ized potential (E2f) becomes complex: 



XV<r,»)=V<p) 5 ; (A.15) 

SV M = - y2(p) ! n(p --> (A.16) 

[i - v(p)»n( P , w )] 2 + [F(p)9n( P , u )f 

Therefore, the functions ^n(p, oS) and ^i2(p, oS) become complex, too, and so does 
the quasiparticle spectrum E{p) = E(p) + ij(p), where 7 is the decrement of damp- 
ing. 

The dynamical structure factor 

S{p,w) = — -3n(p,w) , 

im 

which is measured in neutron scattering experiments and is usually determined 
by the imaginary part of the density- density correlation function n(p, uo) (see, for 
instance, ||43|), can be also obtained as the imaginary part of the "screened" inter- 
action 

s( P|W ) = --^9y(p lW ). 

On the "mass shell", the form factor S(p) = S(p,E(p)) has a maximum in the 
same region of p where the potential V(p, u) and correspondingly the quasiparticle 
spectrum E(p) have a minimum, in accordance with the general conection beetween 



E(p) and S(p) (see Pfl-p 
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Fig. 1. The diagram representation of the Dyson-Belyaev equations for the 
normal Gu and anomalous Gn Green functions of bosons in the superfluid Bose 
system. 
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Fig. 2. The diagram representation of the nonlinear integral equations for the 
normal En and anomalous S12 self-energy parts in the Bose liquid with account for 
terms of first order in the small density of the single-particle BEC (n <C n) and of 
the equation for the renormalized ("screened") potential of pair interaction V. 
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Fig. 3. The radial dependences of the regularized Aziz potential (|39| ) (solid curve) 
and the Lennard- Jones potential fl3?D (dashed curve). Inset, the Aziz potential at 
small distances, down to r = 0. 
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Fig. 4. The finite repulsion potential in the model of "soft spheres" in the form 
of (a) a Fermi "step" ( PH ) (solid curve corresponds to b = 0, dashed curve — to 
b = 0.5a), and (b) a Lindhardt function (|42|) in real space. 
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Fig. 5. The momentum dependences of the Fourier component (pTl| ) of the Fermi- 
type potential (|l0| ) at b = or the Lindhardt-type potential (f42|) (solid curve) and 
of the potential of the "hard spheres" model [^3| V{p) = Vosin(pa)/pa (dashed 
curve). Inset is the Fourier component for the Fermi type potential with b ^ 0. 
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Fig. 6. The Bogolyubov quasiparticle spectrum fl43|) with an oscillating potential 
(p|) (solid curve), maximally close to the experimental spectrum p3|-p^| (dotted 
curve) at Vo/(47ra 3 ) = 169 K and a = 2.44 A. The dashed curve corresponds to 
the unstable spectrum ( K5p at Vo/(47ra 3 ) = 217 K, characterized by negative values 
Eg(p) < in some momentum interval. 
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Fig. 7. The momentum dependence of the renormalized ("screened") potential 
(P^) (solid curve) with account for the momentum dependence both of the Fourier 
component PT| ) (dashed curve) and of the polarization operator U(p,E(p)) on the 
"mass shell" (see Fig. 8). 



27 




Fig. 8. The momentum dependences of the real part of the boson polarization 
operator on the "mass shell", H(p) = 9?n(p, E(p)), multiplied by Vo, obtained from 
self-consistent computations at T = 1 (solid curve), and of V U(p)T at T = 1.5 
(dashed curve). 
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(c) 

Fig. 9. The momentum dependences of the functions (a), ^f(p) (b), A (p) 
(c), obtained from self-consistent computations at the value of the single fitting 
parameter V /{4na 3 ) = 1552 K for a = 2.44 A. 
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Fig. 10. The theoretical quasiparticle spectrum E(p) obtained from the self- 
consistent computations in the model of "soft spheres" for a quantum Bose liquid 
vith a suppressed BEC in the approximation of constant vertex A(p) = T(p, 0) = 1.5 
for all p (solid curve) and the empirical spectrum of elementary exitations in the 
SF Bose liquid 4 He (circles). 
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Fig. 11. The theoretical quasiparticle spectrum E(p) calculated with the vertex 
being a decreasing function falling off from F = 1.5 to T = 1.1 in the region 
2.1 A < p < 3.8 A(solid curve) and the empirical spectra |25 
(asterisks). 
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